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I . Summa ry 

The following report is a description of activities performed and results 
obtained over the past year 12/84-12/85 by The College of William and Mary under 
cooperative agreement (NCCl-91). 

This grant was awarded primarily for developing an adequate technique to reduce 
the Halogen Occultation Experiment (HALOE) data. set. The HALOE experiment will begin 
shortly after launch of the Upper Atmosphere Research Satellite (UARS) in late 1989. 

The successful retrieval of atmospheric parameters from radiometric measurement 
requires not only the ability to do ideal radiometric calculations, but also a 
detailed understanding of instrument characteristics. Therefore a considerable 
amount of time was spent participating in instrument characterization in the form of 
test data analysis and mathematical formulation. 

Both electrical and optical characterization are important for accurate signal 
processing. Analyses of solar-to-reference interference (electrical cross-talk), 
detector non-uniformity, instrument balance error, electronic filter time-constants 
and noise character were conducted. Optical studies included, beam uniformity analy- 
ses, interference between gas-and-vacuum paths (spectral cross talk) and in-orbit 
signal estimates. Consulting roles were played during tests that included field-of- 
view measurements, gas cell response, spectral filter measurements, and signal drift 
measurements to name a few. A second area of effort was the development of 
techniques for the ideal radiometric calculations required for HALOE data reduction. 
The computer code for these calculations must be extremely complex and fast. During 
the year, a scheme for meeting these requirements was defined and the algorithms 
needed for implementation are currently under development. Most of the work in this 
area has been in the development of a new transmittance code. When completed, 

"exact" (commonly referred to as "line-by-line" or "LbL") limb path transmission pro- 
files will be possible in one program execution at previously unattainable speeds. 
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In addition, the formulation for a less accurate, faster radiance calculation method 
has been defined. This method will be combined with the "LBL" method in the overall 
retrieval scheme. 

A third area of work included consulting on the implementation of the EGA (emis- 
sivity growth approximation) method of absorption calculation into a HALOE broad band 
radiometer channel retrieval algorithm. 

Most of these areas of work will be detailed in the following report. 
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II. Introduction to HALOE 

II-l Experiment Overview 

The bulk of this section has been taken directly from the HALOE software release 
I document. It is included here only for the unfamiliar reader. 

The Halogen Occultation Experiment (HALOE) is a satellite solar occultation 
experiment being developed as part of the Upper Atmosphere Research Satellite (UARS) 
program. The objective of UARS is to study the physical processes acting within and 
upon the stratosphere, mesosphere, and lower thermosphere of the Earth. A central 
factor in meeting this objective is the development and launch of a free-flying sci- 
entific observatory to provide needed measurements. Scheduled for launch in 1989, 
this UARS observatory will have a complement of 11 scientific instruments and will 
operate at an altitude of 600 km and an orbital inclination of 57 degrees. As part 
of this instrument complement, HALOE will provide measurements of key species in the 
ClOx, NOx, and HOx chemical cycles (i.e. O 3 . HCl , HF, NO, N 02 » CH^ , and H 2 O), 
and the HALOE team will provide scientific studies with these and related data. 

The HALOE experiment is based on the concept of measuring the absorption of 
solar radiation by atmospheric constituents along a limb-viewing path during sunrise 
and sunset (i.e. occultation) events (Fig. 1). Vertical scans of the upper atmos- 
phere are obtained by tracking the solar disk during an occultation. From these ver- 
tical absorption scans, associated atmospheric temperature data, and orbital informa- 
tion, species concentration profiles as a function of pressure (and altitude) can be 
inferred. This solar occulation technique offers the advantages of high sensitivity 
due to long absorbing paths, high vertical resolution due to weighting of the absorp- 
tion near the ray tangent altitude (Hq in Fig. 1), and a strong relatively constant 
intensity radiation source in the sun. 
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As illustrated in Fig. 2, the HALOE instrument uses gas filter correlation radi- 
ometry for measurements of HCl, HF, CHm and NO and uses broadband filter radioinetry 
for measurement of tangent point pressure (from CO 2 absorption), O 3 , H 2 O, and NO 2 . 

The gas filter correlation concept is based on correlation between the spectral con- 
tent of the solar radiation passing through the atmosphere and that introduced by a 
gas cell containing the species of interest. It provides the high effective spectral 
resolution needed for isolating spectral signatures of a target gas from the effects 
of interfering absorbers. The broadband filter radiometer concept is similar to that 
used successfully on the Nimbus satellites. 

The HALOE instrument provides four types of operations during the science data 
portion of an occultation event: balance, solar scan, calibration, and atmospheric 

data collection. The first three operations occur while HALOE is viewing the sun 
above the atmosphere. In the balance operation, electronic gain adjustments are made 
to match, to some offset, the signals from the two separate paths in each gas corre- 
lation channel. The solar scan mode provides five complete scans (down and back) of 
the sun to characterize the solar relative intensity curve. In the calibration oper- 
ation, a calibration wheel containing gas cells for the four gas filter target gases 
and a series of neutral density filters are cycled through the optical path. This 
provides a check on gas response, radiometric calibration, and instrument balance. 
Finally, in the atmospheric data collection operation, HALOE tracks the sun through 
the atmosphere to obtain modulation (gas correlation channels) and radiance (radi- 
ometer channels ) profiles. For a sunset event, all four operations are included; 
for a sunrise event, the balance is omitted. 

In effect the HALOE technique attempts to differentiate between the absorption 
of the gas of interest and the total broad band absorption, which can be as small as 
one part in 100. This one part in 100 then is desired to be known to at least 10% 
accuracy, requiring knowledge to .1% of the total adsorption. Since the HALOE signal 
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is proportional to transmi ssion (where transmission = 1-absorption) and transmission 
is typically >.9, the total signal must be measured with a precision of .01% or bet- 
ter. As can be seen by the above description, the measurement is done by the differ- 
encing of two large signals. This implies that the consistency of these two large 
signals must be maintained by the instrument and estimated by calculations to better 
than .01%. This accuracy can be obtained only by careful characterization of the 
instrument and high resolution calculations of the spectral correlation between 
atmospheric and gas cell absorption features. 

Both of these general tasks have been addressed over the previous year. 
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III. HALOE Instpjm--^nt Signal 

As discussed above, the success of HALDE will depend on how well the signal can 
be modeled. Therefore considerable effort was given to putting the HALOE signal 
equation into a form that would produce insight and understanding. This would allow 
test results to be interpreted with more confidence, instrument specifications to be 
better evaluated If necessary, and possible ideas to surface on the best attack of 
the retrieval problem. 

HALOE attempts to produce a signal that is highly correlated with the atmos- 
pheric absorption of a specific gas. This is done by effectively differencing two 
occultation radiometer measurements. The input for one optical path measurement has 
spectral radiation, in regions of the specified gas absorption, removed by passing 
the input through a gas cell. The idea is that the two measurements will be similar, 
except for sensitivity to atmospheric absorption by the type of gas in the cell. 

The study of this concept prompted a few questions. 1. Can the signal aV 
(commonly referred to as modulation) be formulated in a way analogous to a radiome- 
ter? 2. If not, what is the difference? 3. Will a change in formulation allow 
radiometer techniques to be applied to the HALOE retrieval problem? 

The HALOE modulation equation is: 

AV = / ^ 

where S = source function (includes normalization (or conversion constant)) 

G = gain factor necessary for balance 
f = broad spectal filter function 
Av = voltage 

Ta = atmospheric transmission 

Tc = instrument response common to both vacuum and gas path 
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Ti = instrument response unique to gas path 
T£ = instrument response unique to vacuum path 
All variables are a function of v (wavenumber) except G which is constant. Also 
V 2 > and f(v) = 0 for v < and v > \> 2 . The integral will be assumed over 
the non-zero range of f(v). 

The instrument is balanced under conditions of = 1, or 



= / Sf(G-Tg-l)dv 

(2) 

where G' 

= Gtj^/t^ 

(3) 

f 

= (Vacuum path instrument spectral response) 

(4) 

One can show that 

(see appendix A) 



AV = AVg - (G - 1) -V^,p + G / Sf'?aegdA 

(5) 


Term 1 Term 2 Term 3 

where V^g = vacuum path signal viewing above the atmosphere 


= vacuum path signal viewing through atmosphere 
^a = 1 - "a 



= lAg 

Tg = / Sf Tgdv/ / SfdV 

Term 1 is set to zero at balance, leaving only terms 2 and 3. Term 2 is propor- 
tional to minus the average broad band absorption while term 3 is proportional to the 
average absorption over the effective gas filter, f'^g (see Appendix A). It is 
interesting to note the broad band absorption term goes to zero for "G = 1. This 
could have been achieved by balancing the instrument with an empty cell in place then 
switching to a filled cell for recording data. This would have allowed the luxury of 
spectral calculations being necessary only where is significantly non-zero. 
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However, the advantage of the HALOE balance is that absorption sources with zero 
correlation to 5g will balance to zero (i.e. wide filter and gas filter average 
absorption are equal). The effects of poorly known continujJ^ or randomly distributed 
spectral features are minimized by the HALOE balance. 

Our first two questions are now answered. Term 2 and term 3 of equation 5 are 
analogous to simple radiometers. As we will see in the next section, term 2 is equi- 
valent to two radiometers with slightly different broad band filters (G'f & f). 

Term 3 of equation 5 is formulated identically to an emission radiometer where the 
filter is G'f^g and the Plank function is replaced by S, the solar source 
function. However, the third question is still in doubt. Although the formulation 
is adequate, the broad band filter calculation techniques may not be capable of 
modeling term 3 adequately. 

To gain more insight into the aV signal, the various components of equation 5 
were evaluated for each of the HALOE gas filter channels using a model atmosphere. 

The calculations were performed for six paths with tangent altitudes of 62, 50, 40, 
30, 20, and 10 kilometers. The computations were simplified by colapsing the optical 
mass into the tangent layer, so the calculations are intended only for a rough study 
of relative magnitudes. The mixing ratio profiles used, are in Table 1. The results 
of the calculations appear in Tables 2a, 2b, 2c and 2d. As can be seen, the terms of 
equation 5 have been normalized by Vyo then divided by 2 x 10"^. This puts them 

5 

into instrument noise level units which are roughly 2 x 10 ” SI^q. The only 
column in the tables that needs explanation is the last. This is an attempt to 
estimate the signal change in Vyi if the principle absorber (the gas of interest) 
were removed from the atmosphere. That is V\l is the evaluation of Vyj 
excluding the gas of interest. Therefore, the last column represents the largest 
possible signal due to the gas of interest. 
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If the interferring absorbers have zero correlation with the gas filter, with 
little overlap of spectral lines, one would expect column 2 and 6 to be nearly 
equal. If there were very little spectral absoption by interferring gases in regions 
of gas filter absorption, then one would expect columns 4 and 6 to be nearly equal. 

By inspecting Table 2a-2d we see that one or both conditions may exist above 30 km, 
but below 30 km there is significant non-zero correlation and/or spectral line 
overlap. This implies that retrieval calculations will need to incorporate detailed 
characterization of line shape at the lower altitudes, especially when Doppler shift 
effects are included. To handle these effects with broad band techniques, large 
empirical adjustments will be necessary and/or a large increase in precalculated data 
tables along with increased complexity in the table utilization algorithms. 


Mi 


ALT (km) HF HCl 

62 2.9E-10 1.3E-9 

50 4.5E-10 1.3E-9 

40 4.0E-10 1.5E-9 

30 3.0E-10 l.OE-9 

20 6.0E-11 5.0E-10 

6.0E-13 2.4E-11 


ng Ratio Profiles Used in S 
TABLE I 

Mixing Ratio Pgag/Pgi 
CHi, NO H2O 

8.0E-9 3.0E-9 3.0E-6 

l.OE-7 l.OE-8 5.0E-6 

4.0E-7 1.2E-8 4.6E-6 

8.0E-7 2.5E-9 4.8E-6 

1.3E-6 2.5E-10 4.4E-6 

1.6E-6 


ations 


N 20 

O 3 

NO 2 

1.7E-11 

0.96E-6 

l.OE-10 

3.1E-6 

3.1 E-6 

l.OE-9 

1.6E-8 

7.3 E-6 

3.5E-9 

9.0E-8 

6.6E-6 

7.5E-9 

2.4E-7 

2.6E-6 

9.0E-10 

3.2E-7 

1.3E-7 

3.0E-10 
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3.0E-10 4.0E-5 


12 


TABLE 2a 
HF Channel 


Tangent 
Ht. km 

AV 

Term 2 
(Eq. 5) 

Term 3 
(Eq. 5) 

V -V 

\o 

vr'^vi 

62 

.11 

-.23 

.34 

2.27 

.16 

50 

.72 

-2.07 

2.79 

21.44 

1.18 

40 

3.39 

-6.27 

9.66 

62.00 

4.52 

30 

11.95 

-16.61 

28.58 

164.10 

13.42 

20 

23.45 

-51.08 

74.53 

504.60 

25.50 

10 

35.23 

-441.69 

476.92 

4362.90 

12.60 

(units of 

2 X 10-^ 

Vvo) = 1* 

1012 




Tangent 


Term 2 

TABLE 2b 
HCl Channel 
Term 3 


'' vr'^vi 

Ht. km 

AV 

(Eq. 5) 

(Eq. 5) 

V -V 

''vo 

62 

.52 

-.02 

.53 

.84 

.52 

50 

2.26 

-.17 

2.43 

9.11 

2.33 

40 

8.74 

-1.64 

10.38 

88.23 

9.43 

30 

28.45 

-12.15 

40.60 

653.80 

33.61 

20 

43.40 

-75.72 

119.13 

4073.72 

81.69 

10 

-45.65 

IT 

-438.55 

393.89 

23592.86 

63.40 


(units of 2 X 10"^ Vvq) ^ 1.0186 
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TABLE 2c 
NO Channel 


Tangent 
Ht. km 

AV 

Term 2 
(Eq. 5} 

Term 3 
(Eq. 5) 

'^vo"''vi 


62 

.05 

-1.72 

1.77 

41.32 

1.41 

50 

16.17 

-7.65 

23.83 

183.00 

21.61 

40 

84.61 

-18.33 

102.94 . 

438.30 

90.78 

30 

164.27 

-49.25 

213.52 

1177.50 

200.00 

20 

94.96 

-176.38 

271.34 

4217.04 

210.76 

10 

42.17 

-948.01 

990.18 

22665.87 

247.66 

(units of 

2 X 10-^ 

Vvo) ^ ~ ^ 

..1418 




Tangent 
Ht. km 

AV 

Term 2 
(Eq. 5) 

TABLE 2d 

CH 4 Channel 

Term 3 
(Eq. 5) 

1 

-V 

'' VI ''v 

62 

.27 

-.25 

.52 

.48 

.22 

50 

14.53 

-10.06 

24.59 

19.09 

17.00 

40 

201.68 

-143.96 

345.65 

273.19 

251.13 

30 

815.02 

-732.17 

1547.18 

1389.27 

1191.21 

20 

2921.96 

-2915.51 

5837.47 

5532.08 

5193.62 

10 

6775.10 

-11749.35 

18524.45 

22293.99 

21596.13 


(units of 2 X 10“^ Vyg) ^ = 1.5270 
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IV. Fast Broad-Band Method - Application of EGA 
(Emissivity Growth Approximation) to HALOE 

Limb geometry remote sensing experiments have the ability to infer unique mixing 
ratio profiles from measured radiance proviles over significant vertical altitude 
ranges, at least to the vertical resolution of the instrument. This vertical resolu- 
tion is determined by vertical fiel d-of-view, sample spacing and noise. For HALOE, 
it should be about 2 km throughout the stratosphere. Because of the extremely non- 
linear nature of the radiative transfer equation, accurate solutions are obtained 
only through iterative procedures involving many successive calculations of the for- 
ward problem. (This inverse problem can be solved in a closed form only if line- 
arized about an atmospheric state assumed "close" to the true solution. This neces- 
sitates an a prior knowledge and/or assumption about the true state, giving unre- 
liable results. (The method is unaccepatable in our opinion.) The ability to do 
many fast and accurate forward calculations becomes the key to credible inversion of 
the measured radiance profiles. 

The most successful method, to date, for doing the broad band forward calcula- 
tion has been the EGA. This method was applied to the LIMS (Limb Infrared Monitor of 
the Stratosphere) experiment operational data processing with excellent results. We 
are now in the process of applying the technique to the HALOE modulation calculation. 

Very simply stated, the EGA method calculates broad-band absorption or emission 
over inhomogeneous paths using tables consisting of broad-band emissivities, 

^(U,T,P), for homogeneous paths. These homogeneous emissivities represent a wide 
range of mass path "U", temperature "T" and pressure "P" values. The broad-band 
emissivity is defined as: 
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r(U,T,P) = / 2F(v)B(v,T)c(U.T.P.v)dv/ f 2F(v)B(v,T)dv (6) 

Vi V2 

where F = Broad-band spectral filter 
B = Plank function 

5 = Monochromatic absorption 

V = wavenumber 

F = 0 outside of range vi to V2 

The use of these tables is described in Gordley and Russell, 1981. 

Equation 6 describes the tables that would be used for an emission calculation. 
For absorption calculations, as for HALOE, the Planck function B(v,T) is replaced 
with the solar source function B(v,6000°K). The same forward emission calculation 
procedures can then be used to calculation absorbed solar radiance. 

From equation A4 of appendix A we see that 

AV = J SfCgClv - / Sf'G'Sadv + / SfG'Cg^adv (7) 

(where AVg is assumed balanced to 0). This is the equivalent of measuring absorp- 
tion using three separate radiometer channels with filters Fj, F2 and F3 where 

Fi = FtcT2 
F2 = GfTcTi 
F3 = GfxcXiCg 

The development of the algorithm to use the EGA techniques in evaluting equa- 
tion 7 is underway with a preliminary functional code likely to be complete by early 


1986. 
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V. Retrieval Scheme 

As already dicussed, the accurate modeling of the HALOE signal must include 
effects that originate from the sharp spectral features in both the gas cell and 
atmospheric spectra. When the relative spectral position or shapes of these features 
change due to Doppler shifting, pressure broadening, etc., the aV signal can incure 
large relative changes. The ability to include these effects in an EGA algorithm 
with consistent and sufficient accuracy over all atmospheric and measurement scenar- 
ios is not only unlikely but difficult to validate. Only exact LbL calculations 
appear to be adequate for this requirement, but their use in interative retrievals is 
prohibitively expensive and time consuming. 

We have developed a scheme that will combine the accuracy of the LbL method and 
the speed of the EGA method while circumventing their respective disadvantages of 
time and inaccuracy. 

Figure 3 displays the retrieval procedure in flow diagram form. The following 
reasoning is the source of this scheme. 

The EGA forward calculations are fast and precise and will be used in the calcu- 
lation intensive iterative retrieval algorithm. However, it is assumed that given 
the correct atmospheric state, the signal levels calculated by the EGA will be in 
error. This error, in our experience, is typically a systematic function of pressure 
and masspath. Therefore, given some measured signal as input, the EGA algorithm will 
relax to an incorrect solution. The LbL algorithm can be used to spot-check the EGA 
error at several target altitudes by using the EGA solution as in^ut to an LbL for- 
ward calculation of signal. The difference between the measured signal and the LbL 
estimated signal is then used to infer a needed adjustment of the signal profile that 
was input to the EGA retrieval program. Stated another way, an adjustment is made to 
the measured signal that effects the correct retrieved solution. This adjustment is 
infered and continually refined by periodic and/or sparse LbL calcuations. 
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Figure 3.- Modulation inversion procedure 
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The measurement adjustment concept may at first seem non-rigorous and unset- 
tling. However, viewed from another perspective, it is in reality in LbL retrieval 
in which an EGA code is used as a very sophisticated relaxation and interpolation 
scheme. To see this, let's look at a likely procedure in determining and maintaining 
an adjustment function "A". 

The problem is to find an adjustment "A" which we'll assume is a function of 
tangent pressure P. 


Me = A(P)M^(P) (8) 

Mni = the measured modulation 

where Me = the modulation that, when input to an EGA code will effect 
the correct retrieval 

Initially we will assume A = 1, retrieve a mixing ratio, then use the mixing 
ratio as input for an LbL calculation of modulation Ml(P). If the EGA calculations 
are reasonably accurate with systematic errors, one might assume that a change of 
AMe(P) will effect a change aMl(P) with the following relationship. 


AMg(P) aM(_(P) 

Using equations 8 and 9 one can show that: 
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A starting scenario may be as follows: 


1. 

Set A(P) 

= 1 and aL(P) = 0 



2. 

Retrieve 

a mixing ratio Q(P) 

using Mg{P) = 

A(P)Mm(P) 

3. 

Stop if 

|A(P)-AL(P) I < e (e 

= convergence 

criterion) 

4. 

Set aL(P) = A(P) 



5. 

Use Q(P) 

for LbL calculation 

of M]_{P) 


6. 

Set A(P) 

= Me(P)/M|_(P) and gc 

1 to step 2. 



The above sequence could involve many LbL executions, depending on relaxation rates. 
However, on subsequent retrievals, the initial adjustment would be done with the 
previously estimated A(P). LbL calculations could be limited by doing only a small 
number of tangent pressures. The update of A(P) at these few pressures could be used 
to infer the other levels through an interpolation or fitting procedure. It may be 
found that even these sparse update calculatons need not be performed every orbit 
since occulation tangent latitudes very slowly with time. There are many possibil- 
ities for holding LbL calculation to a minimum while maintaining retrieval integ- 
rity. The most important aspect of the method is that line-by-line signal simu- 
lations, using retrieval mixing ratios, are required to match measured signals . Th i s 
is the equivalent of a line-by-line retrieval. 
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VI. The LbL Code 

It should be clear from the previous section that the quality of the modulation 
retrieval results will be strongly dependent on the efficiency of the LbL code. 

Speed, flexibility and accuracy are the required ingredients for an adequate LbL 
algorithm. In view of these requi rements , a significant effort was made in designing 
and developing a new Lbl code for this specific application. The goal was to develop 
a code with the following characteristics: 

a) Efficient table look-up Voigt function. 

b) No input/output of cross-sections. 

c) Complete limb signal profile. 

d) Cross-section calculations independent of number of tangent paths. 

Ninety percent of the time and expense of currently available LbL codes can be traced 
to the calculations of the Voigt function and the reading and writing of calculated 
cross section values. Therefore these operations must be accelerated or reduced in 
number before a faster code can be realized. 

VI. - 1 Table Look-Up Voigt Function 

The Voigt function can be viewed as the convolution of two separate symmetric 
functions about a center frequency. These two funtions are a Doppler broadening 
function and a Lorentz (pressure) broadening function. The first is a Gaussion and 
the second takes the form of: 


W[(v-Vq)^ + a^] 
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wnere c is a constant equal to 1 /t, .a is the transition halfwid*'h at STP (standard 
temperature and pressure) and (o - vq) is the wavenumber distance from transition 
canter vg. u, vg and ^ are in units of cm"^. 

The halfwidths of there two functions derive from two different physical phenom- 
ena. The Doppler is due to the distribution of molecular kinetic motion causing 
various doppler shifts of the transition. The Lorentz is due to collisional warping 
of the quantum states causing a change in energy levels. The Doppler broadening is 
therefore temperature dependent while the Lorentz broadening is, to first order, 
density dependent. The range of temperature/density combinations in the atmosphere 
can be extreme, requiring careful planning and validation of any table look-up scheme 
used for the Voigt function evaluation. In addition, the behavior of the Doppler and 
Lorentz broadening functions in the "wings" (large distance in wavenumbers from line 
center) is quite different, exp[-(vg - v)^] versus a l/(vg - v)^ type behavior. 
Therefore, the Lorentz function will always dominate at large distances from line 
center. 

The primary objective was to develop a code that would evaluate the Voigt func- 
tion as quickly as possible with at least the minimum accuracy required. The goal 
for minimum accuracy was set at 1% relative accuracy for points oh the Voigt curve 
with .1% integrated accuracy. These were met with the scheme described below. 

Normally the Voigt function is expressed as a function of two parameters x and 

y- 

X = |v - v^|/ap 
y = a^/ap 

a[_ = Lorentz halfwidth 
ag = Doppler halfwidth 
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If the table look-up is to be fast, the scheme must allow linear interpolation over 
subsets of the tables. The logic to find the proper subset must be fast, implying 
few subsets. It was found that the variation of the function over the x, y plane 
made it difficult to choose table subsets that accommodated easy look-up interpola- 
tions. Therefore a different mapping formulation was used in which the independnt 
parameters were defined as, 

x‘ = |v - VQ|/(aL + ap) 

y' =otL/«o 

Note: y' = y 

x' = x/(y + 1). 

With this transformation of independent variables, it was found that the domain of x' 
and y' could be partitioned into 4 regions, as described in Table 3. 

In the first region a simple Lorentz function calculation is performed. It was 
found that the evaluation of the Lorentz function was faster than table look-up logic 
for extended x' regions. In addition, this allows for distant Lorentz wings. The 
last 3 regions are done using table look-up and linear interpolation between adjacent 
x' and y' values. Figure 4 is a contour of percent error over the domain of x' and 

y'. 


TABLE 3 

TVOIGT Routine Table 
Partitioning 



X' 

Y* 

Limits 

Spacing 

Limits 

Spacing 

Region 1 

X- > Xg 

N/A 

Uni imited 

N/A 

Region 2 

X' < Xg 

.1 

5 < Y' £ 11 

.5 

Region 3 

X 

X 

CD 

.1 

.1 < Y' _< 5 

.1 

Region 4 

I 

CQ 

X 

N/ 

X 

.1 

0 < Y' ^ .1 

.002 


h- 

(11.0 - Y)/(0.6875(Y + 1)) 
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VI. -2 L-b-L Calculation Procedure 

Reducing L-b-L calculation time requires minimizing cross section calculations 
and input/output operations. To achieve this goal with the additional requirement 
that transmission for multiple paths be calculated with a minimum of computational 
overhead, we have ordered and separated the necessary operations in a unique fashion. 

First, the integration over spectral frequency is done as the outside loop. 

That is, all computations involving cross sections over an interval av are 
completed before proceeding to the next spectral interval. This avoids the necessity 
of writing cross-section values to disk or storing large numbers of cross-section 
values for later use. 

Second, the spectral calculations are kept to a minimum by calculating cross- 
sections at pressure and temperature conditions associated with various attitudes of 
the atmosphere being analyzed. It is then assumed that the integrated transmission 
can be modeled sufficiently for each tangent path by a combination of these P, T 
conditions and an associated optical mass. A separate code, given temperature, pres- 
sure and mixing ratio profiles and tangent ray altitudes, will determine the needed 
P, T conditions for cross-section calculations and a corresponding optical mass array 
for each ray path. Therefore, if the L-b-L program recieves N comibinations of (P, 

T) for cross-section calculation and there are M ray paths to be modeled, it will 
recieve an MXN array of optical mass values, for each absorbing gas. 

Figure 5 shows a pseudo FORTRAN code used to descibe the order of calculation. 
Figure 6 graphically depicts the possible relation between levels selected for tem- 
perature/pressure conditions (TPC's) and the mass integration along the tangent ray. 
The masspath is integrated along the ray in finite elements determined by the dashed 
curve borders. The mass of each element is partitioned into the two nearest 
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TPC's. The partitioning is done by a Lorentz half width weighting scheme. For 
instance if T^, and T 2 , P 2 are the temperature and pressures of TPC levels 
1 and 2, then the mass Ue is partitioned as follows; 

Ui = Ug[a2 - ag]/[oi2 - 
U2 = Ug - Uj 

where 


a 


a 


P,- 


1 


= halfwidth at standard pressure and 296 K. 


Y = a gas dependent constant typically .5<y<1. 
Tq,Po = standard temperature and pressure. 


It is obvious that the ray transmission error should approach zero as the TPC's 
become more numerous and closely spaced. This implies that there should be a spacing 
that gives the maximum error that can be tolerated. Studies are planned for deter- 
ming this spacing. Also note that to achieve this desired accuracy at a few selected 
tangent levels, the optimum TPC spacing will likely be nonuniform. 

The code for performing all of these L-b-L calculations is in the final stage of 
development and a prototype should be operational by 3/86. 
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Initialize Arrays 
DO 90 n = 1, N 
DO 50 i =1, I 
DO 40 k = 1, K 

Calculate J cross-section for 

Segment n, SCBi, gask— C(j) 

DO 30 m = 1, L 
DO 20 j = 1, J 

A(j,m) = A(j,l) - u(m,i ,k)*C(j) 
20 CONTINUE 
30 CONTINUE 
40 CONTINUE 
50 CONTINUE 

DO 70 m = 1, L 
DO 60 j = 1, J 

60 T(m) = T(m) + exp(A(j,m)) 

70 CONTINUE 
80 CONTINUE 

N _ _ - jjt of spectral segments 

K # of gases 

L - - - # of raypaths 

I _ _ _ # of spectral calculation boundaries 

J _ _ _ (j! of spectral points in a segment 

C - - - Cross section array 
A - - = Absorption coefficient matrix 
U - - - Masspath matrix 
T - - _ Integrated transmission array 


Figure 5.- Calculation Procedure 



LINE-BY-LINE RADIATIVE TRANSFER GEOMETRY 
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VII. Instrument Characterization Studies 

No instrument is theoretically perfect. Therefore imperfections must be charac- 
terized through testing and modeling. In the past year William and Mary personnel 
have performed studies of optical and electronic "cross-talk", detector non-uni- 
formity, and balance off-set effects. These studies are described in the following 
sections. 

VII-1 "Cross-Talk" Effects 

"Cross-talk" is an expression used to imply the inseparable mixing of signals. 

In the HALOE instrument there are two general mechanisms for creating cross-talk; 
optical reflections and inexact electronic demodulation. 

The first, optical cross-talk is caused by radiation being reflected into 
unintended paths. This causes radiation, nominally intended to fall on a gas 
(vacuum) path detector, to be particually reflected from at least two elements in 
such a way that some significant portion of the reflected energy is directed onto a 
vacuum (gas) path detector. 

The second mechanism, electronic cross-talk, causes detector output due to 
external source radiance (internal black body radiance) to be interpreted as output 
due to internal black body radiance (external source radiance). The internal black 
body is used as a constant radiance source for monitoring the relative system attenu- 
ation between gas and vacuum paths. When the relative attenuation varies, as 
detected by a change in the gas-vac reference signal, a gain is applied to the gas 
path signal by a feed-hack loop called the automatic gain control (AGC). The signal 
due to the internal source is distinguished from external signal sources by chopping 
the two beams of radiance at 300 hz and 150 hz respectively, then electronically 
demodulating the signals. Due to unavoidable imperfections in both the chopper and 
demodulator, there occurs a slight mixing of signals. 
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VII-1.1 Optical Cross-Talk 

Optical cross-talk occurs in a given gas filter correlation channel when a frac- 
tion of the signal from one path (gas filter path or vacuum path) is reflected into 
and interpreted as signal in the other path (vacuum path or gas filter path). The 
amount of error introduced into the modulation signal by optical cross-talk depends- 
upon two factors: the magnitude and the spectral content of the cross-talked signal. 

The HALOE modulation signal, M, is given by: 

M = (GVg-Vy)/VyQ, where 

G = multiplier gain = Vvo/Vgo; 

Vvo = exoatmospheric signal , vacuum path 
Vgo = exoatmospheric signal , gas filter path 
Vv = endoatmospheric signal , vacuum path 
Vg = endoatmospheric signal, gas filter path 

Two types of optical cross-talk were studied. "Gas into vac" (GV) cross-talk 
occurs when signal from a gas filter path is recorded as signal in the vacuum path. 
"Vac into gas" (VG) cross-talk occurs when signal from a vacuum path is recorded as 
signal in the gas filter path. 

If only GV cross-talk is present, (values with a prime denote values with 
cross-talk). 



Vv' = Vv + fVg 

Vvo' " '^VO + foVgo» 

M' = (G'Vg' - Vv')/Vvo' 

fo = fraction of energy on gas path detector crosstalked at balance 
f = fraction of energy on gas path detector crosstalked during event 
Note that f varies with tangent altitude observed. 
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Expressing M' in terms of the V signals without cross-talk, the following 
realtion can be otained: 


M' 



V vn 

VO 0 go 


1 




VO 


+ f V 
0 go 


( 1 ) 


Term 1 is very nearly the modulation, M, with no cross-talk present. Term 2 
represents the source of the greatest modulation error. The fraction of gas path 
energy crosstalked during the event, f, is greater (assuming a non-negative 
modulation) than the fraction of cross-talk during balance, fg. This is because 
spectral features in the atmosphere correlate with specral features in the GV 
cross-talk path. 

The importance of spectral content in GV cross-talk can be seen in the following 
Table 4 profiles simulated for the HCl channel of HALOE. The modulation values are 
expressed in units of noise equivalent modulation, or nem, with 1.0 nem being the 
desired level of agreement. 



TABLE 

4 



GV X Talk 

Effects 


Tan Height 

■ Ml 

M2 

Mg 

62 

.5198 

.520 

.5199 

50 

2.277 

2.278 

2.271 

40 

8.657 

8.658 

8.534 

30 

27.019 

27.022 

26.065 

20 

48.026 

48.045 

42.181 

10 

6.298 

6.362 

-3.623 

Crosswalk level: 

O.Oii; 

. 733 % 

.505% 
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Case has no optical cross-talk and represents the nominal modulation for a 
standard atmosphere. Case M 2 has optical cross-talk of .733% at balance, and case 
M 3 has .505% optical cross-talk at balance. However, the cross-talk path for case 
M 2 was chosen to have a high magnitude but low spectral content. On the other hand, 
case M 3 has a low magnitude but a high spectal content achieved by a cross-talk 
path that makes 3 passes through the methane blocking cell. Clearly, modulations 
M 3 are much differenct from the nominal values M^ than are the values M 2 . 

These cases illustrated the need to know the origin of the cross-talk in the HALOE 
instrument if the modulation is to be corrected properly for cross-talk. 

VG cross-talk can be analyzed in a manner similar to that for gas into vac cross 
talk. For vac into gas, define (again, the prime denotes values with cross-talk) 


V ' = V 

VO VO 


V ' = V + f V 
go go 0 VO 




G' 


M' 


M' 


= V„ 


V + fV 
g V 

V ' 

VO 


go 

G'V ' 

g 


go 


VO 
+ f V 


0 VO 


Vq - 

V„. + f.V. 


which can reduce to 


■go 0 vo VO 

If M denotes modulation without cross-talk, the expression for M' can be 


written 
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M' = M/(l + f,G) +T G(f - fj/(l + f 3) 


■ ? ' 

V '-/ 


1 


where 


V /V 

v' VO 


Here, fg is the fraction of energy on the vacuum path detector crosstalked at 
balance; f is the fraction of energy on the vacuum path detector crosstalked during 
the event. Note that f does not equal fg except at balance. 

Term 1 above is nearly the modulation without cross-talk, as G is on the order 
of .75, for HCl fg is on the order of .002. 

Term 2 indicates the modulation error introduced by the VG cross-talk, and 
predicts that all modulation values with only VG will be larger than the modulation 
values without cross-talk. The Table 5 below lists the results for a VG cross-talk 
case in the HCl channel for .2% cross-talk at balance. 

TABLE 5 

VG Cross-Talk Effects 


t Heights 

Ml (nem) 

(nem) 

62 

0.5198 

0.5198 

50 

2.277 

2.278 

40 

8.657 

8.706 

30 

27.019 

27.423 

20 

48.026 

50.351 


10 


6.298 


9.728 
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Column M| represents the nominal modulations for no cross-talk and column 'Ao 
represents the modulations for .2% vac into gas cross-talk. Note all M 2 entries 
are greater than or equal to M^ entries as predicted by equation 2. Note the 
errors are less than 2 nem above 20 km. 

It should be pointed out that all atmospheric transmission profiles were calcu- 
lated using spectral line parameters from the 1982 Air Force Geophysics laboratory 
tape and line-by-line transmittance code. The HALOE instrument response (i.e., 
optics transmission) was determined by using vendor supplied transmission curves for 
all otpical elements. 

These simulations of the HCl channel demonstrated the need to reduce the magni- 
tude of any GV cross-talk with high spectral content. The HCl channel is unique in 
that it contained a methane gas filter in the optical path immediately following the 
broadband filter. Simulations showed that cross-talked paths making several passes 
through this cell could be responsible for introducing large modulation errors due to 
the high spectral content of these paths. To reduce the magnitude of cross-talk 
paths through the blocking cell, interchanging the location of the blocking cell and 
broadband filter was proposed. Studies of this new configuration showed that modula- 
tion errors due to cross-talk paths through the blocking cell were greatly reduced, 
due to the reduction in cross-talk magnitude achieved by forcing the cross-talk beam 
to more than 2 additional passes through the broadband filter. The configuration 
change was adopted. 

During the course of the cross-talk analyses, it was realized that any cross- 
talk beam will have an effect in both gas and vacuum paths, regardless of its origin. 
This is because any optical train is divided by a beamsplitter into the gas and 
vacuum paths, and there is no direct way for any signal reflecting back out of the 
gas (vac) path to get into the vac (gas) path wihout first encountering the gas-vac 
beamsplitter. Therefore, any and all cross-talk beams are divided into the vacuum 
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and gas paths. This should greatly reduce the error in modulations obtained by 
including cross-talk effects in just on path. In this study we assumed the beam 
splitter would divide cross-talked energy in the same proportion as primary beam. In 
reality the reflected energy will favor the direction from which it came due to the 
polarising effects of the beam splitter. This will further reduced the true cross- 
talk effects relative to our estimates. 

A simulation of gas**vac cross-talk including the component of the beam that 
remains in the gas path was carried out for a high spectral content and high magni- 
tude (4.5%) cross-talk path in the HCl channel. Listed below in Table 6 are the 
modulation values calculated. 

TABLE 6 

Combined Cross-Talk Effect 


Tangent Height (km) 

Mj (nem) 

M 2 (nem) 

M 3 (nem) 

62 

0.5198 

0.5215 

0.5202 

50 

2.277 

2.251 

2.280 

40 

8.657 

7.630 

8.652 

30 

27.019 

18.199 

26.895 

20 

48.026 

-2.869 

47.063 

10 

6.298 

-68.815 

4.100 


Once again, Mi corresponds to the nominal (no cross-talk) modulation level. Case M 2 
is the modulation calculated if only the componet of cross-talk split into the vacuum 
path is included. Note the extremely large changes. Case M 3 is the modulation 
calculated if both components of cross-talk are included. Note that the agreement is 
better than any previous case where only one component of cross-talk was studied, 
especially considering the large magnitude of cross-talk in this last case. 

The entire cross-talk anaylsis was significant indefining the importance of the 
mechanisms (magnitude and spectral content) that contribute to modulation errors of 
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optical cross-talk is present. The study also resulted in the realization that any 
cross-talk beam will add to the gas and vacuum path signals, regardless of its 
origin. Including the cross-talk effects in both gas and vac channels significantly 
reduced the expected error in modulation and also provided a more accurate represen- 
tation of the instuments's expected performance. 

VII-1.2 Electronic Cross-Talk 

The two types of electronic cross-talk are referred to as solar-to-reference and 
reference-to-solar. In general, the solar-to-reference interferes with the A 6 C, 
giving gain correction error and, therefore, aV/V (modulation) error. The reference- 
to-solar causes a constant offset error in the vacuum and gas path signals, also 
resulting in modulation error. For a derivation of the electronic cross-talk 
equations see Appendices B and C. 

VII-1.2.1 Solar-to Reference Cross-Talk 

The error in modulation SgM due to solar-to-reference cross-talk can be 
sufficiently approximated by the folowing: 

53 M = x'/9[e(l - t9) - M*" + mJ] (V-1) 

where x'v vacuum path cross-talk level defined as: 

X V = ^'^vo “ '^rvo^^'^rvo (V-2) 

Note that reference signal voltage is designated by Vp, instead of R. 

^rvo measured vacuum path reference voltage at balance and is the 

m 

measured vacuum path reference voltage while viewing cold space. M is the measured 

modulation, av"’/v"’ . ^ is the measured solar vacuum signal at balance and AV”’ 

svo svo 

is the measured AV signal (in the same units as is measured modulation 

at balance, and can be sufficiently approximated by (aV"’ + 

"e" is a system constant that can be determined by the following equation: 
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X e = 

A. y 


(S - D(1 + »R,/v„„) 


6(H„ * 


1) 


x' M 

^ V 0 


(V-3) 


where 6 is the observed gain, from AGC readout, between balance and space dwell 
(i.e., S = Gj/Gq where = gain at space dwell and G^ = gain at balance). 
is the measured vacuum reference signal during space dwell. aR^ is the balance aR 
value. Again, voltages must be in the same units; i.e., if aR has a gain of 60 
over V|,y, then aR must be divided by 60. Equations V-1 and V-3 can be used in 
combination with balance and space dwell data to completely define a modulation 
correction during track mode due to solar-to-reference cross-talk. 

Preliminary estimates of the effect of solar-to-reference cross-talk can be 
found in Table 7. 

VII. -1-2-2 Reference-to-Solar Cross-Talks 

The error in modulation, 6j^M due to reference-to-solar cross-talk, is suffi- 
approximated by; 


5„H = y,'[(l - t9)cI - if] 


(V-4) 


y^' is the vacuum path cross-talk level defined by: 

. '^vo - V^,(off)/[V^ - V^^(otf)] (V-5) 

svo 

where V*^y(off) is the vacuum path solar measurement while dwelling on cold space 
(off the Sun), and m"’ are defined as before, "d" can be found by 
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d , (V-6 

■^V 

where M"’(off) is the measured modulation while dwelling on cold space. Equations 
V-5 and V-6 can be used in combination with balance and off-Sun data to define a 
correction to the track mode data, necessitated by reference-to-solar cross-talk. 

For a detailed derivation of the electronic cross-talk equations, see appendices 8 
and C. 

VI 1-2. Balance Offset Effects 

As was seen in section III, the HALOE instrument balances the gas and vaccum 
path solar signals so that the voltage difference in outputs is zero. This balance 
operation is performed while viewing unattenuated solar energy, that is, .viewing the 
sun exostmospherically. The signals are matched by applying a gain to the gas path 
signal. After the balance, the solar signals will change due to atmospheric 
absorption during sunset or sunrise. Therefore, in order to maintain this balanced 
condition, an unattenuated internal blackbody reference signal is monitored and 
maintained at its balance level by varying the gain on the gas path signal. This 
reference signal is produced by passing the blackbody emission radiance through the 
same optics as the solar input, to the extent possible. It is chopped at a 
different frequency than the solar input in order to separate solar from reference 
through demodulation electronics. 

The succesful maintenance of balance depends on several factors. First, the 
changes in instrument response to gas and vacuum signals must be the same for both 
reference and solar sources. The variation of system response must not be to fast 
for the AGC to follow. Finally, the required gain changes must not go past system 
capabilities. 
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This last requirement is important when considering initial gain settings and 
balance resolution. Ignoring cross-talk, the reference and solar difference signals 
can be modeled with the equations 


AR = GRg - 


AV = 6Vg - 

where Rg - reference gas path signal 
Ry - reference vac path signal 
Vg - solar gas path signal 
Vy - solar vac path signal 
G - adjustable gain. 

The instrument is initially set, in the lab, so that aR = AV = 0 for G = 1. G 
is allowed to vary, as needed during the mission, from .75 to 1.25. However, aR is 
digitizied and monitored only over a range _+ .03 Ry^ where is the initial 
laboratory setting of R^. 

The balance is maintained over an event by adjusting 6 so that 


G = (aRq + R^)/Rg 

where aR^ is the balance value of aR. If aR^ = 0 it is easy to see that aV^ would 

be maintained at zero for any instrumental signal path response change in the form: 

\ * "g''G 

and/or 


Vy - 
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This is not the case however, if the and/or are different for reference and 
solar signals. One must also consider the balance procedure. At balance, G is 
adjusted to set aV^ = 0, i.e. 


G . (AVj, + V„)/Vg 

Where aVq and Vg are the prebalance exatmospheric values. 

If Kg and vary by 3% or more between reference and solar signals, aR may move 
outside of its digitization limits. The system balancing technique requires that aR 
remain within the digitization limits. One can change these limits and/or balance 
away from 0 (i.e. aV^ * 0) but this will cause necessay corrections to the data. 

One more consideration that should be addressed is what is the impact of an 
imperfect AGO? In other words, can corrections be made for the effects of aR not 
being maintained at the balance condition. 

Appendix D addressed these questions resulting in the following correction 
formula: 


6bM 


6R(M^ + t'') + (Ry"’- Ryjj) (1^- t'' AR'"/r"’^) 
AR^ + if.. 


(V-7) 


V 


Correction to the measured modulation to obtain the true modulation 
measured modulation aV/V, 


VO 


V '"/V 

V ' VO 


VO 

m 


VO 
1 m 


AR' 

5R 


m 


vacuum solar signal at balance 
vacuum solar signal 
vacuum reference signal at balance 
vacuum reference signal 
delta R signal 

(aR'" - aRq) where aR^^ = aR"’ at balance 


Note that equation V-7 contains only measured quantities. 
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Primary points of interest are: 

1. Failure to maintain the balance condition (6R t 0) can be corrected to the 1 NEM 
level since 6R/(R^'^ + aR*'^) should be known to 1 NEM as long as aR stays 
within the digitization range. 

2. The effects of large off-sets in aR at balance are small if 

|(Ry'" - ^v"'H " l*^v “ n small. For example, if aR"^ varied 

by its full range of 0.03, then a 1 NEM or greater effect would require 

|K^ - l| > (2 X 10"^/. 03) - 33 NEM 

The determination of Ky - 1 should be accurate to 5 x 10“** (1/2000) which 
is 25 NEM. Therefore, consideration should be made to extend the range of aR. 

3. Equation V-7 assumes negligible change in the ratio 

^ Kg(ref) K^(solar)^ 

■Kg(solar) Ky(ref) J 

during an event. 

4. If 5R = 0 (A6C maintains aR"' = aR^) then 6gM =0 if = 1. It is insen- 
sitive to Kg. 

VII-3. Correction Formula Interdependence 

The modulation error formulas derived in appendices B, C, and D apply for SRXT 
(Solar to Reference Cross-talk), RSXT (Reference to Solar Cross-talk) and balance 
related errors, respectively. These derivations were done independently; but, in 
reality, are interrelated. The effect of SRXT on the RSXT error formula is 
negligible, as can be seen by formulas C20 and C21. However, all three formulas have 
terms that are dependent on the measured modulation, implying that the corrections 
must be a small fraction of the modulation to be valid when used independently. 
Reviewing the correction formulas, we have: 


SRXT Correction: 
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= x'/^[e(l - t 9) - rf + 


(A) 


RSXT Correction : 

5rM = y‘^[(l - T9)d - M*"] 


(B) 


Balance Related Correction: 


6bM = 


5R(M"’ + t'') + (R^"’ - ■ T''AR"’/Ry"’) 

(aR"’ + R^"’) 


(C) 


An important point to remember is that the SRXT and BR (Balance Related) formu- 
las are caused by AGC driven gain changes (except for the 6R term). One is induced 
by variation in the SRXT signal , the other is induced by system response changes in 
the gas and vacuum paths. However, SRXT changes and system response changes can both 
cause variation in the apparent vacuum reference signal, i.e.. 


K - ‘ » 

This suggests the question, will correction C include any part of correction A? 

The answer is, yes. Since R will change due to the varying SRXT signal during an 

event, R - R ^ will become non-zero without K or K. deviating from 1. This 
v VO 9 V 

should be accounted for in the following way. 
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I'he portion of due to cross-talk is 


x'vRvo(^ - 


To remove the effect of cross-talk from equation C, one should replace it with 
equation DlOA from appendix D 


6 


B 


(aR*^ + Ry"’) 




- 1)] • - t^r/r/) 

(aR"’ + R^*") 


(C‘) 


The validation of these equations must ultimately come from static model calcu- 
lations, lab data analysis, and orbital performance. The first two are in progress. 

One final comment on these correction formulas— it would be aesthetically satis- 
fying if a formula that rigorously couples the three effects could be derived. The 
difficulty with attempts at a general formulation is the complexity of the algebra. 
The cross terms make it difficult to identify terms that can be reliably neglected. 

It is also difficult to put in a form that includes only telemetered parameters. 
Indications from preliminary validation efforts are that a general formulation will 
be unnecessary. 

It should be noted that the time dependence of the AGC response will tremen- 
dously impact the application of the above correction formulas. These formulas are 
valid for direct use only in a static situation, which will be virtually non-existent 
during the mission. However, they serve as tools for understanding the HALOE 
response and basic building blocks for a dynamic correction. 

The point should be made that if all gain changes during an event are determined 
to be false (i.e., due to cross-talk and system balancing imperfections), then the 
aAGC can be used to develop a simple scale factor correction in addition to the 
direct use of equation B. 
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Appendices B, C, and D can also be used to develop more exact corrections, 
should cross-talk and imbalance become large. 

A final point worth considering is a planned calibration sequence for HALOE 
during high S angle (no sunset) times to determine system response changes as a 
function of time and temperature. The full, unattenuated Sun can be observed for 
extended continuous time periods to calibrate system response changes. These 
calibration curves can then be used in combination with cross-talk corrections during 
occultation events. Cross-talk changes as a function of time could be determined by 
periodically scanning off the Sun during these calibration opportunities as well. 

In conclusion, we feel that undesired system response character! sties as a 
function of time and temperature can be determined and monitored in orbit. However, 
it is critical that these responses be well understood, characterized and minimized 
through a thorough instrument test campaign in the lab. 
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Vni. Oetactor to Bean Non-uni forrnity Studies 
Other processes which may induce modulation errors are beam movement on the 
detectors as well as the spatial uniformity of the input beam and detectors. 

Beam movements may induce modulation errors in the following way. If a detector 
has a non-uniform response versus position over its active area, moving the beam to 
another location on the detector may make the signal out of that detector increase or 
decrease, depending on how the average response of the illuminated detector area 
changes. Thus, the same input radiance may give many different outputs if the beam 
moves on a detector having non-uniform spatial response. Beam movements could be 
very critical to HALOE, since the response from two detectors is differenced to 
determine the signal (modulation), and slight (1 percent) changes in responsivity can 
induce a signal error 500 times greater than the noise. 

Detector uniformity data are available in the form of response versus position 
across the active area of each HALOE detector. The data is spaced every 24 microme- 
ters and was produced by a laser spot scan across the detector. The response of a 
detector for a nominal sized input beam (of diameter 360 ym) is just the average of 
the spot scan responses in the area filled by the input beam. 

Studies were done moving the beam a maximum of 24 microns to the left, right, 
above, and below detector center for all four gas filter channels, gas and vacuum 
paths. To study movements of less than 24 microns, the spot scan data were linearly 
interpolated. For simplicity, both gas and vacuum beams were assumed to move in the 
same directions, and by the same amount. A standard radiance and modulation profile 
was used as the nominal response for beams centered on their detectors. To simulate 
the effects of beam movement, the gas and vacuum path radiances were adjusted by the 
percentage response change from the center position to the new location. Modulations 
were then calculated using the new radiances with the multiplier gain calculated for 
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the nominal conditions. Using the nominal multiplier gain assumes the beam movements 
are too fast for the Automatic Gain Control to account for. 

Below, in table 7, is a list of expected modulation values for each HALOE gas 
filter correlation channel at 20-km tangent height, for a movement of 1 micron. The 
center value is the nominal modulation level in NEM. These studies show that small 
(1 micron or less) movements can induce modulation errors on the order of several 
NEM, depending on the direction of beam motion. The methane (CH^)' channel appears to 
be affected least by the movements. This is because the two detectors for this chan- 
nel are more uniform than the other detectors, and both gas and vacuum path detector 
responses change about the same amount for a given movement. 

This study was undertaken in part as an effort to identify possible sources of 
noise observed in the instrument. 

Table 7 

V Signal Change for l-pm Beam Movement 
The following is a list of modulations calculated for a 1-um beam movement: 

91.752 

f 

117.358 + 87.719 ->■ 57.398 NO channel 

81.224 

38.386 

f 

36.327 ^ 48.026 48.498 HCa channel 

46.315 

2988.595 

f 

2987.530 2987.705 > 2987.020 CH^ channel 

2987.425 

34.025 

f 

34.346 ^ 22.981 23.651 HF channel 

+ 

48.169 
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The second study concerning accurate determination of the HALOE modulation 
signals was prompted by laboratory measurements of noise and deals with the possi- 
bility of beam non-uniformities being imaged on the detector. This will occur if the 
instrument is not properly focused. In particular, if beam non-uniformities exist, 
and if those uniformities vary with time, then signal variations may be produced. 

The cause of beam fluctuations would most likely be a fluctuating amount of water 
vapor in the 10-m path between the laboratory radiation source and the instrument. 

The NO channel is most susceptible to water vapor absorption, as a 10-m unpurged path 
absorbs about 14 percent. Absorption in the 10-m unpurged path for the other three 
gas filter correlation channels is about 1 percent or less. Non-uniformities arise 
if the amount of water varies within the beam and across its diameter. It should be 
pointed out that the absolute absorption level is not important in producing signal 
variations but it is the range of absorption over the beam diameter that produces 
signal variations. In addition, the NO channel detectors have the largest spatial 
non-uniformities of the four gas filter correlation channel detector pairs. Varying 
beam uniformity on a non-uniform detector will also yield modulation fluctuations, 
even if the average beam intensity remains constant. 

To simulate beam intensity, non-uniformities and the resulting modulation 
signals, the detector spot scan data was used again. It was assumed that the beam is 
on the center of the detector and that it doesn't move. Instead of taking a 
uniformly weighted average of the detector response over the area filled by the beam, 
a weighted average of the detector response was calculated, with the weights being 
between and 1.0, where T|_ represents the lowest transmittance in the 10-m path 
for a given laboratory relative humidity. The procedure was as follows: 

For each of the four gas filter correlation channels: 

1. Determine for various relative humidities. 
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2. Generate, using a random number generator, a set of weights (beam non- 
uniformities) between and 1.0. 

3. Calculate the weighted average of the gas filter and vacuum path detectors. 

This represents the 'balance* conditions. From these conditions, the ratio of 
the gas and vacuum path responses represents the multiplier gain, G. 

4. Generate another set of weights between Xj^ and 1.0. 

5. Calculate another set of responses for both the gas filter and vacuum path 
detectors. 

6. Multiply the gas filter detector response from step 5 (above) by G and 

subtract the vacuum path response from step 5 (above). This represents the 
quantity GVg - in the numerator of the modulation equation. 

7. Multiply GVg - by 60. (This simulates the magnitude of aV as the HALOE 

instrument would output). 

8. Convert the aV to counts, between 0 and 4095 for aV between ± 5 volts. 

9. Repeat steps 4 through 8 three hundred times to simulate 5 minutes (300 seconds) 
of HALOE telemetry data. 

10. Plot the AV for each simulated 5-minute run. 

Simulations for the NO channel were carried out for Xj^ values of 0.85, 0.9, 
0.95, 0.98, 0.99, and 0.995. The case for x^^ = 0.99 (i.e., an intensity variation 
of 1 percent) closely simulated the amount of noise currently observed in the NO 
channel when the 10-m path is purged with dry nitrogen. The purge path is designed 
to have a maximum 1-percent relative himidity. Note that for the NO channel, a 10-m 
path with one percent relative humidity has one percent absorption. As Xj_ 
decreased, the values of aV increased. 

The CH4 , HCa, and HF channels were all simulated at x^^ = 0.995, since the 
absorption due to water vapor in a 10-m purged path is never greater than about 
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one-half percent. The simulated (SV signals for these three channels produced 
errors within the currently measured noise levels. 

The results of non-uniformity simulations indicate that the NO noise level could 
be produced by a beam with approximately one percent nonuniformities, assuming the 
non-uniformities to be 100 percent imaged on the detector face. However, since the 
aperture stop is intended to be imaged on the detector face, these non-uniformity 
effects are probably upper limits. On the other hand, a less random non-uniformity 
model would produce a larger effect. In addition, the absorption variations over the 
beam are probably overestimated. 

Finally, since lab data taken under various conditions of humidity show little 
change in NO noise levels, beam non-uniformity variation is probably not the source 
of the noise observed in the NO channel. 
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IX. Consultation Services 

The principal and associate investigators have also served as consultants in 
many other areas of the HALOE project. Both algorithm development and test review 
and analysis were involved. A brief description of these tasks follows: 

IX-1. HALOE Instrument Testing Support 

The past year has seen a comprehensive battery of characterization tests 
performed on the HALOE instrument. Interpretation of the eventual data set will be 
critically dependent on the instrument characterization determined by these tests. 
Therefore, it was necessary for the investigators under this grant to be closely 
associated with the evaluation of test specifications, procedures, and results. 

Specifications were evaluated or derived by estimating the signal sensitivity to 
the particular instrument function being characterized and requiring that the know- 
ledge of that function allow signal modeling to be accurate to 2 parts in 10^ 

(1 NEM). 

These evaluations were done using the HALOE spectral instrument model and/or error 
studies like those presented in appendices B and C. 

Test readiness meetings were attended as consultants for help in determining if 
the test procedures would meet characterization requirements. 

Finally, test results were analyzed to determine if objectives were met, both in 
instrument characterization and performance. Of considerable utility in this process 
was a plotting code developed by the associate investigator on this grant, for the 
display and statistical analysis of the HALOE telemetry stream. 

Some of the tests involving William and Mary personnel participation were: 

1. Filed-of-View Characterization 

2. Empty Cell Bias 

3. AV Bias 

4. Balance Linearity 
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5. Cross-Talk (electrical and optical) 

6. Field-of-View Matching 

7. Dynamic Range Setting 

8. Calibration Wheel Characterization 

9. Command Mode 

10. Throughput Measurements 

11. Gas Response Test 

IX-2 Algorithm Development Consultation 

The successful processing of the HALOE modulation channel data will require 
careful attention to the preinverson processing of the telemetered data. The removal 
of instrument effects, pressure registration, and solar liml) darkening curve measure- 
ment must be done precisely and accurately. The same is required of the radiometer 
H 2 O channel retrieval since H 2 O profiles are necessary for interferrent correction in 
the modulation channels. The techniques and/or code being implemented to perform 
these functions originate from work performed by the principal investigator of this 
grant, although predating this grant. Therefore, the principal investigator's assis- 
tance to individuals implementing these functions for HALOE has and will continue to 
be a significant portion of the contract support. 
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X. Conclusions 

The past year has seen significant progress toward the ultimate goal of success- 
fully inverting the future HALOE data set to desired mixing ratio profiles. 

The primary accomplishment has been the detailed definition of the modulation 
retrieval approach. Besides the definition, critical software modules required by 
the approach are nearing completion, namely, the Line-by-line code and associated 
modules for fast Voight function evaluation, interpolation, path integration, and 
line parameter logistics. 

In addition, the application of the EGA technique to production signal simula- 
tion has been outlined in detail. The software needed for implementing EGA is in 
hand. The implementation procedure has been broken into steps equivalent to broad 
band uses allowing completion of the task in about 2 months. A prototype modulation 
retrieval package should be operational by early spring of 1987. 

A second area of accomplishment was the mathematical characterization of the 
major signal response errors. Balance offset, aR drift, optical cross-talk, solar- 
to-reference cross-talk, and reference-to-sol ar cross-talk have all been quantified 
with in-flight correction procedures defined. It should be noted, validation by the 
instrument static model and lab data analysis are still needed and will be performed 
in the coming year. 

The final major area of support was in consultation on instrument testing. This 
involved many hours of meetings as well as test result analyses. 

The coming year should see the successful software completion, testing, and 
documentation of the modulation channel retrieval algorithm. 

Consultation support in the areas of instrument testing and HALOE system soft- 
ware will continue as before. 



Appendix A 

Modulation Formulations 


IV = / SfTjt^CGTgTj - Tj)dv 

f = fx^Tj 
G‘ = Gtj/t2 

Equation A1 becomes 

AV = / Sf Tg{G'Tg - l)dv 

Balance equation is 

AVg = J Sf (G'xg - l)dv 
^a “ ^ ~ ^a 

Equation A2 becomes 

AV = / Sf (1 - Ca)[G'(l - 5g) - l]dv 


or 

AV = / Sf [G'(l - Cg) - l]dv + / Sf Cgdv - / Sf SgG'dv + / Sf 

Si nee 

/ Sf [G'(l - 5g) - l]dv = AVg 


we have 

AV = AVg + / Sf - I Sf CgG'dv + / Sf SgCgG'dv 
for illustrative purposes note that if 


G' = constant = G 
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which is likely a good approximation, then 


/ Sf^gClv - / Sf'CgG'dv = (G - 1) / Sf (1 - Tg)dv 

Defining 

VyQ = Vacuum path measured signal exoatmospherical ly 

Vyj = Vacuum path measured signal viewing through atmosphere 

Then equation A4 becomes 


4V = AVg - (G - 1)(V^Q - Vyj) + G / Sf ^a^gdv 


Modulation is normally defined for the HALOE experiment as: 

M = AV/Vyo 

From the balance equation we see that for AVg = 0 (as with HALOE) 


IT = 


Sf'dv 


TF 




= i/7„ 


Equation A5 can be seen to be 


(A5) 


(A6) 


where 


M = 


AV 

^0 


-(- 


1 - X, 


4(1 - 




(1 - 




(1 - 




(A7) 


Tg = average wide band transmission, or as weighted by Sf'//Sf'dv 

= average gas filter transmission as weighted by Sf'(l - x )//Sf'(l - x )dv 

g g 

This points out that the HALOE signal is proportional to the difference between the 
average absorption over the gas path filter defined by f'^^ and the broad filter 
defined by f. 



Definition of terms in this appendix not previously defined 
f = broad band spectral filter 
= atmospheric transmission 

= attenuation common to both gas and vacuum paths 
Ti = attenuation unique to gas path 
T2 = attenuation unique to vacuum path 
S = source function 

G = gain factor applied to gas path to meet balance condition 



Appendix 3 

Sol ar-to-Reference Cross-Talk 




The electronic processing of the HALOE signal from a detector can be modeled as 
follows (assuming system linearity): 

V = ojRFxQ (Bl) 

where 

V = voltage 
u = watts on detector 

R = responsitivity of the detector in amps/watt 
F = conversion factor of volts/amp 

X = efficiency at which the a.c. signal is converted (demodulated) to a d.c. signal 
Q = a fixed gain, assumed constant 

For simplicity we can combine Q and F, assuming they are fixed. 


V = (jRxC (82) 

The reference blackbody power on a gas channel detector, interpreted as reference 
signal, would give a voltage of: 

'^rvg “ *^rg^rg^rrg^vg (33) 

where subscript r indicates reference, rr indicates reference input to reference 
signal and g indicates gas channel. Likewise, the reference signal due to solar 
input in the gas channel would be: 


'^srg “ ‘^sg^sg^srg^'rg 


(B4) 
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Let 


and 


Then 


rg ^rrg'^rg 
■^g ~ ^srg'^^rrg 



0) 

u 


D 

rg^rg 

sg^sg 


(B5a) 

(B5b) 


(B5c) 



V + V 
rrg srg 


= A C + Y A C * 
rg rg sg rg 


(B6) 


In reality HALOE does not make a measurement, therefore, this term will even- 

tually be eliminated from our formulas. 

The analogous equation for the measured reference signal in the vacuum path 
would be: 


where 


/"* = 
rv 


V +V =AC +yAC 
rrv srv rv rv '^v sv rv 


Xy “ ^srv^^rrv 


(B7) 


(For the ideal instrument, x^rg = Xgpv " 

For this development we will assume no reference-to-sol ar cross-talk; i.e., 
Xpjy = 0. This gives the solar signal equations of 


Xrsg 


V 

V 


m _ 

V . = 

A C 


sg ^ 

ssg 

sg 

sg 

m _ 

V... = 

A C 


sv 

ssv 

sv 

sv 


(B8a) 

(B8b) 


We now define fs and fp as: 

*Superscript m indicates a measured quantity 



58 


f. = 


'"sg 


f_ = 


c 

rg 

rv 


al so 


It can then be shown that 


V = A C = V 

rg - ^g^'rg Vrg 

v = A r' = V 

rv ■ rv rv rrv 


ar"’ = - v"!^ 

rg rv 


r ' f 
rv r^g r,.,m 


V .V + V 7^ [^(V" - V*" )] 

Vg rv ^v TT^y'-Xyl^s'' sg sv^-i 


In order to simplify the algebra we will make the following definitions; 


(B9a) 

(B9b) 


(BlOa) 

(BlOb) 


(BID 


Xv = ^v^'rv^^sv 
I i (Xgf,)/(X/S) 

Designating a balance condition by subscript o equation Bll at 

< = Vgo - Vvo * X,(IV^g„ - 

At some later time the measured solar signals change while aR^ 
applying a gain we will call G to the gas channel signal 


(B12a) 

(B12b) 

balance would be: 

(B13) 

is held constant by 


ar"’ = aRq 



rvo 


X (IV"" 
^v' sg 



(B14) 
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Note that is a constant voltage while is a changing measurement that 

includes both G and atmospheric effects; i.e., Subtracting equation 

B14 from B13 and solving for G gives: 


[V 


rgo 


X (V^ 
^v' sv 


) 

svo' 


X I(V^ 

^v ' sgo 


- v"’ )]/V 
sg^-*' rgo 


(B15) 


The error in modulation 6M is the difference between measured modulation less modu- 
lation for G = 1. Remember, G is implicitly contained in v'T ♦ 


6M = (V^ - v"’ l/V*^ 

^ sg sv^^ svo 


(V”g/G - = (GV^g - V^g)/(GV^„„) (B 16 ) 


Substituting B15 into B16 gives 


5M = 


.,m " 
^sg^v 


GV'" V 
svo rgo 


[^V*" - v"’ ) 

^ sgo sg^ 


(v'" - v"' )] 

' svo sv'-^ 


(B17a) 


Defining 


e = I - 1 


(B17b) 


and assuming M*'’ = at balance, equation B17 becomes 


A 

To determine the magnitude of this error, x^> e and G must be measured or 
calcul ated. 

In the lab, a cross-talk value Xy is measured by balancing the instrument 
while viewing a laboratory black body and recording a measurement for The 

input is then shuttered-off causing the which is recorded. Xy is then 


defined as: 
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x' = fv"’ - V )/V 

' rvo rvo^ rvo 


from equation B13, this can be seen to be 


(B19) 


x' = V v"’ /V 

svo' rvo 


(B20) 


Substituting B20 into B18 gives 


rgo 

^ sgo 


x'V 


m 


sgo sg'' sgo 


(B21) 


Defining = V^g/V^go and setting = 1, equation B21 becomes 

6M c. x;T9[e(l - t 9) - Nf" + mJ] (B23) 


Two tasks remain before B23 can be used. The validity of B22 must be estab- 
lished (this will be done later) and e must be determined for each channel, e can 
be found directly if one scans off the Sun after balance, dwelling on deep space 
until the AGO (automatic gain control) drives aR back to aR^ by applying an 

A 

additional gain of G. From equation B13 we would have, off the Sun: 


aR"* = GV - = aR„ 

rgo rvo o 


(B24) 


Multiplying equation B13 by G, substracting it from B24, and solving for I gives 


_ ^svo (G - 1) ^^rvo ^^0^ 

m ^ m ^ 


sgo 


V*" X 

sgo^v 


(B25) 


Since 
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and from B20 


B25 goes to; 


= 1/(1 + f^) = 1 - 
svo' sgo O' 0 


X = x'V /v'" 

A. Ay pyQ/ 5yQ 


I - 1 = e = 


(G - 1)(1 H- 4R„/V,,„) 
Gx;(l + M^) 


(626) 


From the AGC telemetry one can determine G while equation B19 can be used to deter- 

/S Ak 

mine Xy* Note (G - 1)/G = AG incurred from on to off the Sun. 


One should note from equations 812a and B20 that 

x' = V C /(C V ) 

Ay Ay'-py 'svo' ' SV rvo^ 


(B27) 


This shows that the cross-talk value Xy “ V^yo (other parameters stay nearly con- 
stant regardless of input). This means that the orbital cross-talk effect will be 
larger than experienced in the lab. Also note that 


^ Xg^p XgCpgCgy 


(B28) 


>.y s 


or from B8a through BlOb 


I = 


V V A A 
Ag’prgo rvo sgo 

” yV A 

rrvo rgo svo 


(B29) 


From B29 we see that I is constant, so we will find a relation for a perfectly 
balanced system (i.e., = 0 and aR^ = 0). One can show from equations B13 and 


P^- 1 - (x;,e) 


rvo 


B19 that: 
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for Xy « 1 and Xg = Xy . then 


A A 
rvo sgo 

A A 
rgo svo 


(B30) 


when the system is balanced. The A's are amps out of the detectors at balance due 
to reference and solar watts on the detectors. These ratios can be obtained directly 
from Bob Spiers throughput measurements and are included in Table Bl. Note that the 
assumption of Xg = Xy is not reliable and, therefore, I should be determined by 
dwelling off the Sun. 

TABLE Bl.- CROSS-TALK PARAMETER ESTIMATES 


Channel 

^svo'^^sgo 

A /A 
rvo' rgo 

e 

^v(2680) 


^v(6000) 

NO 

.5650 

.5714 

.011 

.0003 

3.41 

.001 

CH^ 

1.272 

1.2187 

-.042 

.0024 

4.13 


HCt 

.7695’ 

.7485 

-.027 

(-.040) 

.0028 

4.17 

mm 

HF 

.7855 

.7676 

-.023 

.9943 

5.31 

B 


The Xy(2680) measured in the lab. The Xy(5ooo) value is the 

predicted level of cross-talk that will be encountered in orbit. It is also esti- 
mated that e for the HCt channel will change from lab estimates due to the differ- 
ence in the solar and lab Planck curves as weighted over the gas cell. The solar 
value is enclosed in parentheses. Other channels are predicted to incur little 
change. 

One final task is to show the validity of B22. This is done by first observing 


that equation B14 is 
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aR*" = GV + Y (IV^ _ ) - V 

rgo sgo sv^ rvo 


(B31) 


or from B20 


GV 


. (■^s.o - . .R' 

V™ „ ''rvo 

svo 




rgo 


,m 


rvo 


1 - 


(B32) 


rvo rvo 

Using values in Table 2, Table Bl, and a balance within aR range of 1/60 of 


rvo 


^^rgo , , ,,__g , aR"’ 


- 1 < x'(ex9 - Nf") + All— < .016 


rvo 


rvo 


This proves that 


GV 


rgo 


rvo 


1 


V 

rvo 

m 


1 


rvo 


and that the error in 6M due to this assumption is less than 1.6 percent of 6M. 

Since this appendix was written, it has been determined that Xy can be quite 
different for the aR demodulation. This negates the use of B19 for the determina- 
tion of Xy» However, XyS can be determined by the use of B26 since from B26 


x;(e + M^) x;e = 


G . 1 (1 > ARo/Vpvo) 
: (1 + Mj 


(B33) 


For a balance near zero, the correction equation B23 can be approximated by 


5M = XyST^Cl - T^) 


(B34) 
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In addition, if equation B30 is used in conjunction with throughput measurement of 
e, B33 can be used to infer x'y* However, it is hoped that B34 will be an adequate 
correction, which will be the case if cross-talk does not increase dramatically from 
lab levels. 

For the static case, these corrections have been verified with the static model, 
but under the assumption that the aR and R cross-talk mechanisms are identical. 
The dynamic case will cause a greater effect. Since the gain change, when input is 
blocked, is 6 - 1 = Xy6» the error in gain at sunrise will be x'e* For low beta 
angles, the complete event may be only one or two AGC time constants long, preventing 
the gain from returning to balance conditions. From Table Bl, we see that this can 
be on the order of a 20 NEM error, from lab measurements. However, the Xy values 
are currently larger than listed in Table Bl due to the aR versus R demodulation 
differences as mentioned above. It is hoped that setting the demodulation phase by 
minimizing aA 6C change, when input is shuttered, will improve the effective cross- 
talk. In any case, a time dependent correction will be needed. This may be as 
simple as a direct use of the aAGC, if other gain change causes are negligible. 



Appendix C 

Reference-to-Sol ar Cross-Talk 
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In order to follow this analysis, a thorough review of Appendix B is required. 
Equations analogous to B3 through B56 would be: 


M - W R X r 
ssg " sg'^sg'^ssg^sg 

'^rsg '^rg'^rg^rsg^sg 

*^sg ■ ^ssv^sg 

^g ^ ^rsg/^ssg 

Y = X /X 
V rsv' ssv 

The equations analogous to B6 and B7 would be: 


= 

V.. 

+ V 

A C 

+ Y A C 

sg 

ssg 

rsg 

sg sg 

g rg sg 

= 

V 

+ V 

A C 

+ Y A C 

sv 

ssv 

rsv 

sv sv 

V rv sv 


(Cl) 

(C3) 

(C4a) 

(C4b) 


(C5) 

(C6) 


Equation Bll and B12 equivalents would be: 


c:.. Y„f, 




y 

V 


= y C /C 
V sv sv 


(C7) 

(C8) 


Defining 


Equation C7 becomes 


J = (Yg1^s)/(Yv^r) 


AV"* = 


sg 


- V 


sv 


y (JV 
J'v' rgo 


- V ) 
rvo^ 


(C9) 
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If J is very nearly 1 and d 5 J - 1, at balance 


AV'" = V - V + y [(1 + d)V - ] 

0 sgo svo ‘ rgo rvo-* 


Now define H such that 


A\/" = HV - V 
0 sgo svo 


Subtracting Cll from CIO and solving for H 


H = 1 + Z/V 


sgo 


where 


Z = y [(1 + d)V - V ] 
V rgo rvo 


The measured modulation without cross-talk would be 


M = [(1 + Z/VggQ)Vgg - VgyD/VgyQ 


Measured modulation with cross-talk: 


= (V - V + z)/v"’ 

' sg sv ''svo 


The lab measured cross-talk y^ is; 


y' = (v"’ 

Jy ' svo 


V )/V 
svo" svo 


C15 can be written as: 


m"* = (V - V + Z)/V - y'm"’ 
' sg sv " svo •'v 


The error in modulation is (C17 - C14) : 


SM = If - M = Z(1 - V^g/V3g„)/V^,„ - y;M” 


(CIO) 


(Cll) 


(C12) 


(C13) 


(C14) 


(C15) 


(C16) 


(C17) 


(CIS) 
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Defining C18 becomes 


- n V - V + dV 

6M = y (1 - r9)f rgo rvo rgo ) _ 


svo 


Since ^19 becomes 


M = y'[(l -xSk^IAvo*^) . rf>] 


rvo 


rvo 


From equations B13 and B20, one can show that 


V - V 
rgo rvo 

V 

rvo 


= v' fr 
A. \/ 


svo 


V m 
svo 


X 'e 

A-v 


and, therefore 


V 


- -'e + 1 


= Xv 


rvo 


From Table Bl, it is apparent that C20 can be approximated by 


6M = y;C(l - x9)d - m""] 


Observe from definitions of J and I that if Xg = Xy 5 yg = J “ 1 
then d = -e. 

Table Cl shows y^ values for each channel, both lab and estimated 
levels. The d values are simply the negative e values from Table Bl 


(C19a) 

(C20) 

(C21) 


(C22) 
and I == 1, 

orbital 
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TABLE Cl.- REFERENCE ^ SOLAR CROSS-TALK 


Channel 

y; 2680 K 

y; 6000 K 

d 6000 K 

NO 

-.003 

-.0088 

-.011 

CH 4 

.002 

.0005 

.042 

HCii 

.002 

.0003 

.040 

HF 

.0025 

.0005 

.023 


Note that for NO is negative. This is not unexpected. Demodulation phase 
errors can cause negative signals. Also note that y^ is inversely proportional to 
input radiance levels. 

In order to find d operationally, use the following development. At balance: 


,m 




AV = V‘“ - V“' = V - V + y (JV - V ) 

0 sgo svo sgo svo •'v' rgo rvo' 


(C23) 


Off the Sun: 


A A 


AV*" = l/” - = y (JGV - V ) 

sg sv "^v' rgo rvo' 

ar"’ = = GV - V = aR 

rg rv rgo rvo o 


(C24) 

(C25) 


where is the same gain described in Appendix B, that which must be applied to 
hold AR = aRq. 

From the definition of y^ 


y; = v^^( 0 ff Sun)/(v;\^ - V^^(off Sun)) 


(C26) 
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from C25 


A 

V = (aR + V )/G 
rgo ^ 0 rvo''^ 


Substituting C27 into C24: 


iv"(off Sun) = y,[J(AR„ + 


From C19b and C28, making use of C19b 


AVm _ Vq) 

^ v'^svo ^rvo 


(C28) 


(C30) 


Now making use of the relationships 

V = /(y‘ + 1) 

svo svo'^ 

and 

V = /(x‘ + 1) 

rvo rvo ''^v ' 

We have, using the definition of J = d + 1 

N"(y; + 1) ARo<*v * 



(C31a) 


(C31b) 


(C32) 
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of if the balance is such that aR^ is within the digitization range of 1/60 (v'J1^q, 
"d" can be approximated to 1.6 percent by: 


for expected 
estimates of 


vro 

levels of aRq, Xy. and y^; i.e., ARo/'^rvo ^ ^ ^ 
cross-tal k. 


(C33) 

and current 


. „ [^(off Sun) 

(J s . 

•^v 

or 

y^d 2 !^(off Sun) 


(C34) 


(C35) 


A simpler way, going directly from C30 and C16 (since 




(off) = = 


m 
'svo 


)P” 


+ 1 


[J(# 


+ 1 ) - 1 ] 


rvo 


Since d = J-l,.then 

lT(off) =-^ ( ^ + d) =. y;d 
jy ^ rvo 

We should point out that the simplest way of correcting for this reference-to-sol ar 
cross-talk error would be to subtract the M"'(off Sun) value from the measured 
modulation. Any resulting balance offset is treated accordingly. The above 
development is useful in determining the magnitude of the expected modulation error 
if cross-talk is neglected. 
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At balance 


APPENDIX 0 

aR Offset Effect on Modulation 



(Dl) 


where 

aRq = measured gas-vacuum reference signal at balance 
R = measured vacuum path reference signal 

VO 

Rgo = gas path reference signal 
At some time after balance 

aR^ ^ aR*" = GK„R„„ - K„R„„ (D2) 

0 g go V VO ' ' 

where 

Kg = a system response change factor in the gas path 
Ky = a system response change factor in the vacuum path 

G = a gain applied to the gas path signal in an attempt to maintain aR*" at 
balance value. 

Using Dl and D2 to solve for K^: 

Kg = (aR'" + R;J)/CG(aRq + (D3) 

where 


Ry “ *^v*^vo measured quantity) (D3a) 

If Kg = K^ = 1, then no system response change occurs and the true modulation would 
be measured as: 

= (Vg - V,)/v;^ 


(D4) 



= solar gas path signal 
= solar vacuum path signal 
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V ^ = V = measured vacuum path signal at balance. 

VO VO r 3 

Assuming system changes occur, and that they have the identical multiplicative effect 
on solar and reference signals , then the measured modulation would be (Note that 
cross-talk is ignored.): 


= (®gVg - 


(05) 


We would like to determine a correction factor 6gM = - M^, in terms of measured 

quantities, in order to correct the measured modulation to true modulation values. 


6 n = = If - 

''vo 


(D6) 


Using D5 to solve for and substituting D3 for K_ . we have 

g - a 

''g = (""''vo * * ’'vo'''!''"'" 


(07) 


where 


Vv = 


(D8) 


Now D6 becomes (using D7 and D8) 


b';(m”Vvo v‘;)(aRq ^ r^q) - v';R',o(ar'" I- R';) 
R'Jv^JiR"’ + r'J) 


,fm< 


•mr 


6bM 


= M*" - 


(D9) 


Continuing the algebra with the substitution of aR^ = aR*^ - 6R (which defines 6R) 


j „ = r';[>^Vvo(ar” - .SR) ^ ^ v!;(ar^ - 5R) ^ v!;r,„ - v!;r,„] - v;R',„&r'" 

® * 'C> 
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M^6R 


- \o> - - "»") 


VO' 


RyV^6R 


(ar"’ + r'J) 


r"’v (ar"’ + r"’) 

V VO ' V ' 


R>vo("«"’ " ^v) 


Further reduction along with the identify 


t'' = 


v VO 


gives 


R*^ - R 

= (fT + t'') + — — (m"’ - t''aR"’/r'") (DIO) 

(aR"* + R'J) aR*" + Ry '' 


Since r'" - R is partially due to solar-to-reference cross-talk it must be 

V VO 

accounted for. The portion of r'J that is due to cross-talk is 


The portion of R^ - R^^ that is due to cross-talk is 

- 1) 


and DIO becomes 


6 bM 


6R 


(AR"* + r'J) 


(M^ + t'') + 


K - ''vo 


(DlOa) 


AGC error 


aR' + R^ 

Error due to aR balance offset 


Note that the correction accuracy limitation for balance offset comes from 

one-half bit of r’’’ and R . If R*” is digitized 0 to 2000 then for 1 NEM 

V vo V 

accuracy (We have ignored component and denominator error.) 
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2x1/2 V aR'" 
2000 


2 X 10“^ 


If = 0.5, we must have 



0.08 


before a 1 NEM error effect is incurred. 

It is useful to show what happens if system response changes do not afftect both 
solar and reference signals identically. To do this, we assume that there occurs a 
system response that affects only the reference signal; i.e., 02 becomes 


aR"’’ = GK'R„„ - K'R„„ 
g go V VO 


(DU) 


05 would be 


= (GVg - V,)/V„ 

In this case and = v'JJ 

6M' = 1^' - = m'"' - (Vg - V^)/VyQ 


( 012 ) 


(013) 


Using the same algebra as before: 

6M' = (G - 1)(m"’' + t'') (014) 

In this case G should be monitored using the AGC data. From Oil and the 
balance equation, 01, it can be shown that 


G = 




K ‘(Aft + 

g ' 0 




(D15) 
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Or, if aR"’ = aRq = 0 



(D16) 


G is the ratio of the changes in system response for vacuum and gas reference 
signal, but only when the solar signal responses are unchanged. 

A more general description (although not yet adequately developed) is as 
follows. Let the various response change factors be: 

Kgg = solar gas 

K „ = solar vacuum 

vs 

K „ = reference gas 
gr 

K = reference vacuum 
vr 

Then 


s ■ Ss 

(017a) 

K, > K,s 

(D17b) 

K = K„K„' 


gr g g 





or 

Kg' . Kg^/Kgj (D18a) 

K/ = K„p/K„ (D18b) 

If all four individual response change factors were known, then D14 could be used to 

calculate a new M*" = m"’' - 6M' using G obtained from D15 or 016. A new = 

R'"'/K ' could be calculated. Then the new aR"™, and r") could be used in DIO 
V g V 

for a further correction to the modulation - 6gM. However, it is unlikely 

that these four response factors could be quantified. Besides, the only likely 
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mechanism for causing Kg' abd K^' to be other than unity is independent solar 
solar and reference beam movement on the detectors during an event. In addition, if 
Kv' = Kg', then from D15 we see that G would be unity and D14 would give 6M' = 0. 
Therefore, equation DIO should be sufficient for corrections due to instrument 
response changes. 

One final note. It could also be theorized that response changes may occur for 
only the solar signal. Since this would exhibit no change in AGO, it would be diffi- 
cult, if not impossible, to detect— much less quantify. However, equations D18, 016, 
and 014 could be used to make the correction if the response factors were known. Of 
course, in this case, equation 010 would not apply since no gain change occurs. If 
these changes were determined to be a function of time from acquisition, correction 
equations might be deduced from analyzing the data of long duration exoatmospheric 
tracking events. This idea will be pursued. 


